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We study the mixing structure of isospin-singlet scalars, the light quarkonium (qq) and glueball 
{gg) in two-flavor QCD, based on a holographic model of bottom-up hard-wall type. In the model 
the pure quarkonium and glueball states are unambiguously defined in terms of the different U{1)a 
charges in the restoration limit of the chiral U{2)l x D(2)_r symmetry, in which the quarkonium 
gets massless as the chiral partner of the pion. Hence the qq-gg mixing arises in the presence of 
the nonzero chiral condensate or pion decay constant. At the realistic point where the pion decay 
constant and other hadron masses reach the observed amount, we predict the tiny mixing between 
the lightest quarkonia and glueball: The smallness of the mixing is understood by the slightly small 
ratio of the chiral and gluon condensate scales. The low-lying two scalar masses are calculated to 
be ~ 1.25 GeV and ~ 1.77 GeV, which are compared with masses of /o(1370) and /o(1710). Our 
result implies that /o(1710) predominantly consists of glueball. 


I. INTRODUCTION 

The spectrum structure of the low-lying scalars around 
the mass 1 GeV still remains unsolved in the low-energy 
QCD, and has currently been pursued extensively by sev¬ 
eral approaches. Of particular interest are the isospin- 
singlet scalars, denoted as /o’s in the quark model [I|, 
since they can arise as mixtures having the same quan¬ 
tum number = O’*"*', such as the light two-flavor 

{u,d) quarkonium [qq), glueball [gg). Understanding 
such a rich low-lying isoscalar spectrum is therefore of 
great importance, which is tied to revealing some part of 
the dynamics of low-energy QCD. 

The straightforward investigation for the mixing struc¬ 
ture has so far been performed on the full lattice QCD 
simulations [^, which focuses on the mixing qq and gg 
stateSjand also has been made by using QCD sum rules 
(Ref. [3| for a recent review). The quantitative estimate 
on the mixing angle is, however, still less accurate so 
that one cannot say anything about the constituent struc¬ 
ture. On the one hand, other recent approaches based on 
phenomenological hadron models have attempted to give 
some insight on the qq-gg mixing. By performing global 
fit to model parameters with use of the phenomenological 
inputs such as currently observed scalar-decay properties, 
it has been indicated that /o(1710) is almost constructed 
from the pure glueball state jj. Still, however, the issue 
is controversial so it needs more indication from some 
different approaches. 

In this paper, we study the mixing between the quarko¬ 
nium andglueball states based on the gauge-gravity du¬ 
ality d, dl, so-called holographic (^D. We employ a 
holographic model proposed in Ref. [7| , an improved ver¬ 
sion of bottom-up hard-wall type proposed in Refs. ii- 
The model has succeeded in reproducing the character¬ 
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istic features of QCD, such as the observed meson in¬ 
teraction properties as well as the ultraviolet-asymptotic 
behavior of QCD Green functions. What should be re¬ 
marked in the model of [3| is, in particular, simultaneous 
incorporation of the five-dimensional bulk fields dual to 
the gluon and the chiral condensate operators, and 
qq. As emphasized in [ 3 ], one can thereby reproduce the 
ultraviolet scaling of QCD current correlators including 
terms along with the condensate of G^^ as well as the 
leading logarithmic term and the chiral condensate term. 
The isospin-singlet scalars, arising as fluctuation modes 
around the condensates in the holographic bulk, thus in¬ 
clude both the quarkonium and glueball states, hence the 
mixing among them can be evaluated straightforwardly. 

We analyze the two-flavor QCD case, in which the chi¬ 
ral U[2)l X U[2)ii symmetry is spontaneously broken 
down to the vectorial one U[2)v via the nonzero vacuum 
expectation value of a bulk-scalar dual to the qq, accord¬ 
ingly due to the nonzero pion decay constant /^. In the 
model the pure quarkonium and glueball states are unam¬ 
biguously defined in terms of the different U[1)a charges 
in the restoration limit of the chiral U[2)l x U[2)n sym¬ 
metry [{qq) and /^ —>■ 0), in which the quarkonium gets 
massless, reflecting the chiral partner of the pion. Once 
the chiral condensate develops from zero, the qq-gg mix¬ 
ing is turned on and grows monotonically with the chiral 
condensate. At the realistic point where /,r and masses 
of other mesons such as the vector/axialvector mesons 
reach the desired amount, we find the tiny mixing be¬ 
tween the lightest quarkonium and glueball, which is con¬ 
sistent with the large Nc picture on that of basis the 
holographic model has been established. This is the def¬ 
inite prediction obtained without any phenomenological 
inputs such as observed isospin-singlet scalar decay prop¬ 
erties, in contrast to other effective-hadron model ap¬ 
proaches We further provide a new insight on 

the small mixing in terms of generic quantities in QCD: 
the smallness of the mixing is understood by the some¬ 
what small ratio of the chiral and gluon condensate scales 
at the realistic point. 



2 


At the realistic point the lowest-lying two scalar masses 
are calculated to be ~ 1.25 GeV and ~ 1.77 GeV, which 
are compared with masses of /o(1370) and /o(1710). The 
lightest isospin-singlet scalar is almost completely degen¬ 
erate with the isospin-triplet scalar, ag meson, which 
is due to the extremely-small mass split induced from 
the vanishingly small mixing among the isospin-singlet 
scalars. These two lower masses ~ 1.2 GeV can actually 
be lifted up to the desired value around ~ 1.3 — 1.4 GeV, 
when the mi xing with a four-quark state would be taken 
into account | 1 ^ . 

Our result thus implies that /o(1710) predominantly 
consists of glueball with the mass around ~ 1.7 — 1.8 
GeV, in accord with the pure glueball mass estimate by 
lattice simulations [ 3 , a recent study based on a holo¬ 
graphic QGD of top-down type and a different ap¬ 
proach based on a phenomenological model J. 

This paper is organized as follows: In Sec. n]we give a 
brief review of the holographic model in Ref. 7| and sum¬ 
marize things necessary for analysis of the isospin-singlet 
scalar mixing. In Sec. Illll we define the pure quarkonium 
and glueball states and introduce effective interactions 
describing the mixing of them deduced from the present 
holographic model. The effective mass matrix, obtained 
by keeping a few low-lying scalars, is then evaluated in 
Sec. IIVI Summary of this paper is given in Sec. [V] 


II. MODEL: PRELIMINARIES 

We begin by reviewing the holographic model proposed 
in Ref. [Tf and list some equations and formulas necessary 
for the later discussions. 

The model in Ref. Q is based on deformations of a 
bottom-up approach for successful holographic dual of 
QCD [ 1 , The model we shall employ is described as 
U{2)lX C/(2)h gauge theory which is defined on the five¬ 
dimensional anti-de-Sitter (AdS) space-time. The five¬ 
dimensional space-time is characterized by the metric 
ds^ = gMNdx^dx^ = {L/z)^ {rj^^dx^dx’' — dz"^) with 
= diag[l, —1, —1, —1]. Here, M and N {g and i/) 
represent hve-dimensional (four-dimensional) Lorentz in¬ 
dices, and L denotes the curvature radius of the AdS 
background. The fifth direction, denoted as z, is com- 
pactified on an interval extended from the ultraviolet 
(UV) brane located at z = e to the infrared (IR) brane 
at z = Zm, i.e., e < z < z^. The UV cutoff e will be 
taken to be 0 after all calculations are done. 

Besides the bulk left- (Am) and right- (Rm) gauge 
fields, we introduce two bulk scalars and $gg. The 
4>gg transforms as a bifundamental representation field 
under the U{2)l x U{2)r gauge symmetry, and there¬ 
fore is dual to the quark bilinear operator qq having the 
U{1)a charge equal to 2. The ^gg is, on the other hand, 
chiral and C/(l)A-singlet, dual to the gluon condensate 
operator The mass-parameter for these two bulk 

scalars and are then holographically given as 

Ml = —‘ilL^ and Ml = 0, which reflect the scaling 

SPqg ' ^gg 


dimensions for qq and respectively. 

The action of the model is thus written as Q 

S5 = ‘S'bulk + S'uv + 'S'lR , ( 1 ) 


where 


^buik = [ d'^x dz^^ 

J Je 9b 


UM'i'gg'X 'V 


gg 


+Tr[GM4>JgG^cI>gg - 


QQ 
■MN 


+RmnR^^^] 


( 2 ) 


-S'lR = d*x dzS{z - Zjn)\/^Cm, 


( 3 ) 


with the boundary-induced metric g^i, = {L/zYr]g,v 
The covariant derivative acting on <I>gg in Eq.([2]) is de¬ 
fined as DM^qq = djA^qq + iLlA^qq “ i^qqRjA, where 
Lm{Rm) = Ll,{R%)T'^ with = fTVy2 (a = 0 , 1 , 2,3) 
being the generators of SU{2) and ( 7 ° = 12 x 2 normalized 
as Tr[r“T^] = <5“^ L{R)mn is the five-dimensional field 
strength which is defined as L{R)mn = dML{R)N — 
dNL{R)M — i[L{R)M 1 L{R)n], and g is defined as g = 
det[( 7 MAr] = (A/z)^°. The UV boundary action S'uv in 
Eq.d!]) plays a role of the UV regulator to absorb the UV- 
divergent e terms arising from the five-dimensional bulk 
dynamics, which we will not specify. The IR boundary 
action A'ir is introduced so as to realize minimization of 
the bulk potential by nonzero chiral condensate Q with 
the IR Lagrangian: 

Air = (-m^Tril^ggp] + ATril^ggfl^) . ( 4 ) 


The gauge coupling 35 and a parameter c appearing in 
the action are fixed by matching with the UV asymptotic 
behavior of vector/axialvector current correlators as Q 


L _ _ L _ 

g\ 127r^ ’ 1927r3 ’ 


where the latter has been determined from terms damp¬ 
ing along with the gluon condensate. 


A. Vacuum expectation values of scalars 


The bulk scalar fields $gg and are parametrized as 
4>gg(a;, z) = ^ {vqq{z) + (Jgg{x, z) -b al{x, z)ct*) 

(i = l,2,3) (6) 


X(a:, z) = e = v^{z) 


( 7 ) 


with the vacuum expectation values, Vgg = \/ 2 ($gg) and 

cg| . . 

= ix) = e 2 We hereafter disregard pion fields 
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7r“ which will not be relevant for the present study. Plug¬ 
ging these expansion forms into the bulk action 5'buik 
Eq.Q we find the coupled equations of motion for the 
vacuum expectation values Vqq and v^, 




1 


,3 ‘'99 


vld^Vqq I —Vqq = 0 


dz ( -^dzVx I + 2 s L 


z 

2 r2 


{dzVqqf 



where 



1 


- 0 . 01 . 


0 ( 8 ) 


(9) 


The last equality of the second line in Eq. © follows from 
Eq.®. The boundary conditions for Vqq and are cho¬ 
sen (in the limit where the sources for qq and are 
turned off) as 0] 

^99 (^) “ 0 , Vqqi^Zm) — ~j~ j 

%(e) = l) v^{zm) = l + G, (10) 


with the IR boundary values ^ and G which are holo¬ 
graphically related to the chiral and gluon condensates 
as shown in Ref. Q . Looking at the numerical value of s 
in Eq.®, as done in Ref. Q one may neglect the second 
term of 0 {s^) = G(10“^) in the equation of motion for 
in Eq.® to decouple the Vqq and Then the ap¬ 
proximate equations of motion can be analytically solved 
to be Q 


V, 


99 


^t(1 + G) 


{z/Zmf 


L'-^ • l + G{z/zmV’ 
+ G{zlZm)'^ , 


( 11 ) 


with the boundary conditions in Ea. ifTOll incorporated. 
[We have checked that, by numerically solving the cou¬ 
pled equation in Eq.® keeping the term, this approx¬ 
imation is satisfied within a few percent level in a wide 
range of the parameter space we will consider later. ] 

As noted above, the nonzero value of the chiral conden¬ 
sate parameter ^ can be realized by adjusting the IR po¬ 
tential in S'lR of Eq. O: putting the solutions in Eq. m 
into the action S 5 in Eq. © and performing the mini¬ 
mization of the potential energy, one finds the stationary 
condition Q 


f = 


mf,L — 


Nc 

127r2 


3- 


4G 
1 + G 


( 12 ) 


B. Scalar sector 

The five-dimensional bulk scalar fields aqq{x,z), 
cFgg{x, z) and ao(x, z) in Eqs.® and ® can be expanded 
in terms of normalizable modes, corresponding to QCD 
scalar mesons: 

(Jqq{x, z) = Y^ a''^\x)fi^2^z ), 


^99 ( a ^> ^) = Y1 ’ 

n 

a}){x,z) (13) 

n 

with the wavefunctions , fi^l , /i"^ having the nor¬ 

malizable UV boundary conditions, 

/W(e) = 0, /i:He) = 0, /^:](e) = 0. (14) 

Using Eg. (ITS)) we expand the bulk scalar sector of the 
action Sbuik in Eq. © to find the induced effective La- 
grangian in four-dimension: 

^scalar = + £^2 -|- £^^2^ -|- , (15) 


where 



In deriving the effective Lagrangian we have used the 
equation of motion for Vqq and in Eq.® and imposed 
the IR boundary conditions for and /a"] so as 
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to eliminate the IR boundary terms in quadratic order of 
fields as follows: 

fi:l{Zm)=0, (18) 

where use has been made of the stationary condition for 
Vqq in Eg. (1121) . It should be noted that the qq-gg mixing 

strength in Eq.((l7l) is proportional to the chi¬ 

ral condensate ^ Vqq. This implies that the qq-gg mix¬ 
ing is turned off when the chiral symmetry is restored, 
where one can realize the definitely pure qq and gg-isospin 
singlet-scalar states, as will be discussed in the next sec¬ 
tion. 


247r^ 

247r2 


Ar + 3 - 


-AC^ + S- 


4G 

TTg 

4G 

TTg 




C. Vector and axialvector sectors 


The five-dimensional vector and axial-vector gauge 
fields Vm and Am are defined as 


Vm 


Lm + Rm 


Am = 


Lm — Rm 

V2 


(19) 


It is convenient to work with the gauge-fixing 14 = 
Az = 0 and take the boundary conditions I4i(x, e) = 
Vfiix), Afj,{x,e) = a^{x) and d^V^{x, = 

dzAfiix, z)\z=Zm = Oj where v^{x) and afi{x) correspond 
to sources for the vector and axial-vector currents, re¬ 
spectively. We then solve the equations of motion for 
(the transversely polarized components of) V^(x,z) and 
Afj^(x, z) and substitute the solutions back into the action 
in Eq.d!]), to obtain the generating functional IT[n^,a^] 
holographically dual to QCD. Then we obtain the vec¬ 
tor and the axial-vector current correlators, which are 
defined as 


i / dTe*^"(0|TJ“;(4(x) J^^(0)|0) 


= (5 

with the currents 


ab ( 




( 20 ) 


\ 


= g- -= 7T 


V2j 


JT = qi A 


V2j 


7 759- 


( 21 ) 


Av{Q^) and XVa{Q^) (where Q = \/—q^ is the Euclidean 
momentum) are expressed as 


Av[Q^) = 


2 , A dzV{Q^z) 


127r2 


2=£—>0 


AA{Qn = 


2^ _ Nc dzA{Q^,z) 


127r2 


( 22 ) 


where the vector and axial-vector profile functions 
V{Q‘^,z) and A{Q^,z) are defined as Vf^{q,z) = 
Viz{q)V{q'^) and A^^z) = a;i(g)A(g^) with the Fourier 
transforms of v^{x) and a^(x). These profile functions 
satisfy the following equations: 

[-Q"^+ w-^dzu:dz\V{Q'^,z)=Q, (23) 


-Q^ -I- uj ^dzUjdz - 2 



A{Q^z) 


0, (24) 


w = -n" , (25) 

with the boundary conditions V{Q'^,z)\z=e^o = 

A(Q'^,z)\z=e^o = 1 and dzV{Q'^,z)\z=z^ = 

dzA{Q^,z)\z=Zm = 0 - 

The vector and axial-vector current correlators, Tly 
and IIa, can be expanded in terms of towers of the vec¬ 
tor and axial-vector resonances. We then identify poles 
for IIv.^ as the infinite towers of p and oi mesons. Their 
masses, mp„ and are calculated by solving the 

eigenvalue equations for the vector and axial-vector pro¬ 
file functions Q: 

Vn{z) = 0, (26) 


’^(ai)n + ^ ^dzUjdz - 2 



Ajii^z^ 


0, (27) 


with the same boundary conditions 14 (e) = 4.„(e) = 
0 and dzVn{z)\z=z^ = dzAn{z)\z=z^ = 0 . Note that 
the Pn and (ai)n meson masses get degenerate when the 
chiral condensate ~ Vqq ~ ^ is sent to zero. 

The pion decay constant is expressed in terms of 
Ily and IIa as = IlyfO) — n^(0). We can express this 
/. by using Eq.dSH) as 0: 


A dzA{Q\z) 

127r2 z 


2 =e-s -0 


(28) 


The vector and axialvector sectors are completely fixed 
once the parameters 4 G and Zm are chosen to be certain 
values. In Refs. 0 and the optimal values are found 
so as to reproduce the experimental values of Z^, masses 
of the lowest vector and axialvector mesons rup^ = rrip 
and = ixiai- the optimal numbers leading to the 

realistic point are 


Z~3.1, G~0.25, zT~ 347MeV. (29) 


This parameter choice yields /,r — 92 MeV, rup ~ 775 
MeV, TOqj ~ 1264 MeV 0 and predicts 


i-qq) ~ (277MeV)3 , cz (331 MeV)^ , 

TT ^ 


in agreement with values estimated in Refs. 



^ Our estimate on the chiral and gluon condensates is also consis¬ 
tent with values obtained based on other QCD calculations [l3l 


z—e^O 
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III. REALIZATION OF PURE qq AND gg 
STATES AND THEIR MIXING 

As noted above, in the present model the qq-gg mixing 
term is proportional to the chiral condensate ^ Vqq [See 
Ea. (fT7)l ]. When one considers the ideal limit where the 
chiral symmetry is restored {{qq) —t 0, /^ —>■ 0), there¬ 
fore, the pure qq and gg-sta,tes can be unambiguously 
defined. In the present model this limit can actually be 
achieved by sending the chiral-condensate parameter ^ to 
zero [See Ea. dTT]) ]. 


A. Meson mass spectra in the restoration limit of 
the chiral symmetry 


In the chiral-restoration limit, the masses of the pure 
qq and gg isospin-singlet scalars are determined by the 
following decoupled eigenvalue equations derived from 

-Cscaiar in Eq.((Tni) with ^ = 0: 






3 


1 + 




Ksfo)) 


[& = 0, (31) 


where the superscript 0 denotes ^ = 0 {{qq) = fn = 0)- 
Note that, in the chiral-restoration limit the wavefunc- 
tion for the isospin-triplet scalar Oq"^ obeys the same 
mass eigenvalue equation as that of the isospin-singlet 
scalar (Jqq'^ ■ Since the same boundary condition is im¬ 
posed on [fa^l]° and as in Eas. lfHl) and (fTSl) . the 

isospin-triplet and -singlet quarkonia are completely de¬ 
generate in the limit ^ 0 {{qq) —0 and /,r —t 0). 

In addition, the vector and axialvector mesons satisfy 
the same eigenvalue equations in the limit ^ ^ 0 [See 

Eas. (l26l) and (l27l) ]. so their masses are degenerate as well. 
Note the explicit independence of ^ on the eigenvalue 
equation for in Eg. (1^51) . We may therefore fix the 
parameters G and Zm to the optimal values in Eg. (1291) so 
that the p meson mass already reaches the realistic value 
in the chiral-restoration limit ^ 0, i.e., 

m° = = mp ~ 775 MeV, (32) 

where, again, the superscript 0 stands for ^ = 0. 

Taking ^ = 0, G ~ 0.25 and ~ 347 MeV, we thus 
calculate the scalar meson masses in the chiral restoration 
limit. [One should note that the dependence of the po¬ 
tential parameter A in the IR boundary condition Eo. dTSl) 
goes away when ^ = 0.] The lowest-three scalar masses 
below 2 GeV are found to be 


~ 1782MeV, 

~ 1782MeV. (33) 


The lightest quarkonia, for both isospin-singlet and - 
triplet scalars, thus become massless in the chiral- 
restoration limit, reflecting the chiral partner of pions. 


B. Defining mixing between pure qq and gg states 


Once the chiral condensate develops from zero, the 
isospin-singlet qq and gg states start to mix according 
to the mixing form given in Ca a in Eq. (1161) . After the 
scalar fields are canonically normalized, the scalar mass 
terms in Eg. dTSl) take the form 


C 


scalar 




2 1 
1 
2 


k^n 


where 


TDTl 

V C7gg / 7 

TDn 

V'^ao ) Jin ’ 

flo 
on 

^ Ggg) J^n ’ 




g; 


{k,n) 


Ai JA 


(35) 


Note, however, that the scalar fields and are no 
longer purely qq and gg-staies due to the nonzero mixing 
triggered by the nonzero 

For the purpose of addressing the mixing between the 
pure qq and gg-states, we shall now propose an effec¬ 
tive mixing term inspired by the present holographic 

(n) 

QCD. We first replace the mixed wavefunctions /A,, and 
fa^l in the off-diagonal mass-squared element 
of Ea. (l35|) with the pure qq and gg wavefunctions in the 
chiral-restoration limit (^ —0), and [/a"]]® in 

the decoupled equations in (I5T]) . Then, the off-diagonal 
mass-squared element (mt^g"^)^ is modified as 




^{k,n) 


Afe ./A™ 




= 2 V 2 . 






(36) 


where Vqq = L ■ Vqq and (A) = J dz{L/z)^v‘^{z)A{z) for 
an arbitrary function A(z). Thus the redefined (mi^’g"^)^ 
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as in Eq. (1361) is evaluated as the mixing amplitude of the 
pure qq state overlapped with the pure gg state, which 
are properly defined in the chiral-restoration limit, to be 
developed by nonzero ^ only through the nonzero Vqq. 

On the other hand, we keep the diagonal elements 
and unchanged and assume that the 

wavefunctions in these elements satisfy the decoupled 
equations similar to Eg. (1311) in the chiral-restoration 
limit, but with the different IR boundary condition for 
aqq'^ having nonzero ^ as in Ea. dTS)) . 

We have also computed the mass eigenvalues directly 
by solving the mixed equations, and checked that these 
drastic deformations do not make significant difference in 
the scalar mass spectrum. This implies that the mixed 
wavefunctions and fa^l are well saturated by the 
pure wavefunctions [fa^]]° and in the scalar mass 

spectrum. 

Even when the chiral condensate grows from zero, the 
isospin-triplet scalars Oq are still isolated as well as in 
the case of the chiral-restoration limit, except the dif¬ 
ferent IR boundary condition with nonzero term in 
Ea. (fT51) . We may fix the IR potential parameter A so that 
the lowest oq mass is set to the desired value at the real¬ 
istic point achieved by the parameter choice in Eg. (0^1) . 
As done in Ref. , for the reference value of the oq mass 
we may choose mag — 1.24 GeV by taking into account 
that the mass can be lifted up to that of ao(1450) when 
the mixing with a four-quark state is incorporated M- 
Then, the optimal value of A is fixed to be Q 

A = 7 —^ • K , with /c ~ 1. (37) 

(47r)2 

IV. ANALYSIS OF THE MIXING STRUCTURE 

We are now ready to discuss the mixing strength be¬ 
tween the pure qq and gg states, which evolves from zero 
at the chiral-restoration limit = 0) to the realistic 
point ~ 3.1) with the parameter choice in Eas. (l2^ 
and (1571) . In Eig. [1] we first show the evolution of the di¬ 
agonal mass-squared elements and {m'aggY 

respect to /tt, in place of Here, we have only picked up 
the low-lying three mass values, 'm'a}g{= mag), m^gg and 

ml^gg . The figure tells us that the lowest-qg scalar mass 
dramatically develops from zero at the chiral restoration 
point {fgr = 0) to the fixed mag — 1-24 GeV at the re¬ 
alistic point (/^ ~ 92MeV). The growth of the mass 
will actually be saturated around the value ~ 1.3 GeV, 
even if one extremely takes the large fg^. The lowest-gf; 
scalar mass migg is completely independent of /^, while 
the mass of the first-excited gg-state, m^Jg , goes beyond 
that of the ground state of the gg-scalar, to be above 2 
GeV at the realistic point, but not to be above 3 GeV 
even when /^ —>■ oo limit. 

In Fig.[2]we display the mixing strengths (jn^gp)'^ and 


im^gg'^)'^ as a function of fgg. We see from the figure that 
the mixing strengths are still small enough, even after 
reaching the realistic point, compared to the diagonal 
mass-squared values: the mass-squared matrix in unit of 
GeV^ looks like 


/ (1.25)2 (0.14)2 Q \ 

m2 ~ (0.14)2 (177)2 _(o.12)2 , (38) 

\ 0 -( 0 . 12)2 ( 2 . 26)2 J 

which acts on the vector (crqj\ , Cgg^)^. By diagonal¬ 
izing this matrix, one obtains the masses of the low-lying 
two isospin-singlet scalars below 2 GeV, 

mi ~ 1.25 GeV, m 2 ~ 1.77 GeV. (39) 


The lowest mass can be lifted up to the amount of 
the /q(1370) mass ~ 1.3 — 1.4 GeV in a way similar 
to the an (1450) case through mixing with a four-quark 
state [^. Thus, compared to the diagonal element in 
Eo. (l38l) . we find that the mixing strengths are negligibly 
small, in accord with the assumption made in Ref. [7| 
and the large Ng picture on which the holographic model 
has been based. 

The smallness of the mixing can be understood by 
somewhat smaller chiral-condensate scale than the gluon 
condensate scale: these condensates are actually related 
involving the parameters G and Zm as Q 


32 G {-qq)jzm 

3-\/3 1 + G (^Gfgg,) 


(40) 


Taking G ~ 0.25, {^Gl„) ~ (331 MeV)^ and z-^ ~ 347 
MeV, we numerically evaluate ^ to write 


^ ~ 5 X 


(-gg) 

(331 MeV)3 


~ 5 X 



(41) 


with the chiral and gluon condensate scales Agg = 
(—gg)^/^ and Agg = At the realistic point, 

from Eo. (|3nil we have (Agg/Agg)^ ~ (0.84)^ ~ 0.6, in 
accordance with the optimal value of f in Eg. (1551) . If 
(Agg/Agg)^ wus US mucE as/more than of G(l), ^ would 
get larger to enhance the mixing strengths to be compa¬ 
rable with the diagonal element in Ea. (l38ll . as expected 
from Eig. [2] 

Our result thus implies that /o(1710) is predominantly 
constructed from the glueball state with the mass around 
~ 1.7—1.8 GeV, in agreement with the pure glueball mass 
estimate by lattice simulations 0, a recent study based 
on a holographic QGD of top-down type 0 and a dif¬ 
ferent approach based on a phenomenological model Q . 


V. SUMMARY AND DISCUSSION 

In this paper, we studied the mixing structure of 
isospin-singlet scalars, the light quarkonium (gg) and 
glueball (gg) in two-flavor QGD, based on a holographic 
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m5r*_j*^(solid),m®^(dashed) and (dotted) 



FIG. 1. The diagonal mass-squared elements 'rn^^^{= mao) 
(solid curve), (dotted curve) and (dashed curve) 

versus the pion decay constant with the parameter G, 
and A hxed to be the optimal values, G — 0.25, = 347 

MeV and A = 3/(47r)^ (See Eas. (l2^ and (l37t ). 



/dMeV] 

FIG. 2. The same for the mixing strengths (mig’g^)^ (solid 
curve) and [—(mi^g’g^)^] (solid curve) as Fig. [T] 


model of bottom-up hard-wall type. Based on that 
model, we have demonstrated that the qq-gg mixing takes 
place due to the nonzero chiral condensate. In the model, 
the pure quarkonium and glueball states are unambigu¬ 
ously defined in terms of the different U{1)a charges in 
the restoration limit of the chiral U{2 )lxU (2)/{ symme¬ 
try, in which the quarkonium gets massless as the chiral 


[1] S. Eidelman et al. [Particle Data Group Collaboration], 
Phys. Lett. B 592, 1 (2004). 

[2] A. Hart et al. [UKQCD Collaboration], Phys. Rev. D 74, 
114504 (2006) hep-lat/0608026 . 

[3] W. Ochs, J. Phys. G 40, 043001 (2013) arXiv:1301.5183 
[hep-ph]]. 


partner of the pion. We showed that, at the realistic 
point where the pion decay constant and other hadron 
masses reach the observed amount, the lightest quarko- 
nia and glueball are hardly mixed at all, without any 
phenomenological inputs such as the currently observed 
isospin-singlet scalar-decay properties. This is consistent 
with the large Nc picture on that of basis the holographic 
model has been established. The smallness of the mix¬ 
ing strength can actually be understood by the slightly 
smaller ratio of the chiral and gluon condensates (See 
Eg. (ITT]) ). The low-lying two scalar masses are calculated 
to be ~ 1.25 GeV and ~ 1.77 GeV, which are compared 
with masses of /o(1370) and /o(1710). Our result implies 
that /o(1710) predominantly consists of glueball. 

Several comments are in order: 

In addressing the scalar meson masses, we have incor¬ 
porated not only the lowest masses from each of qq and 
gg states, but also the next-to-lowest one from the qq 
state. As seen from the mass matrix element in Ea. (l38l) . 
the mixing between the two lightest scalars from the qq 
and gg could be comparable with that between the light¬ 
est scalar from the gg and the next-to-lightest one from 
the qq. This would imply some new possibility that in¬ 
corporating the next-to-lowest qq scalar could be relevant 
in investigating the low-lying scalar mixing structure. 

Also, we have neglected the higher level more than 
the third in the quarkonium state and the second in the 
glueball state. Actually, we have computed the higher 
masses to be > 3 GeV, which has not yet been established 
in experiments on the scalar resonances, but could be 
accessible in the future, together with the slightly lighter 
qq-like scalar with the mass ~ 2 GeV predicted in the 
present model (See Eg. (1551) ). 

Similar analyses as done in the present paper can be 
performed in a more realistic situation, where mixings 
with a four-quark state and a pure strange state with the 
sizable strange quark mass are incorporated by extending 
to the three-flavor case, which will be pursued in another 
publication. 
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